New insights on the matter-gravity coupling paradigm 
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The coupling between matter and gravity in General Relativity is given by a proportionality 
relation between the stress tensor and the geometry. This is an oriented assumption driven by 
the fact that both the stress tensor and the Einstein tensor are divergenceless. However, General 
Relativity is in essence a nonlinear theory, so there is no obvious reason why the coupling to matter 
should be linear. On another hand, modified theories of gravity usually affect the vacuum dynamics, 
yet keep the coupling to matter linear. In this Letter we address the implications of consistent 
nonlinear gravity /matter coupling. The Eddington inspired Born-Infeld theory recently introduced 
by Banados and Ferreira provides an enlightening realization of such coupling modifications. We 
find that this theory coupled to a perfect fluid reduces to General Relativity coupled to a nonlinearly 
modified perfect fluid, leading to an ambiguity between modified coupling and modified equation of 
state. We discuss observational consequences of this degeneracy and argue that such a completion of 
General Relativity is viable from both an experimental and theoretical point of view through energy 
conditions, consistency, and singularity-avoidance perspectives. We use these results to discuss the 
impact of changing the coupling paradigm. 
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Introduction. - Within Einstein's theory of General 
Relativity (GR), the existence of singularities, e.g., at 
the big bang or within black holes, seems unavoidable 
[1] . Worldlines of freely falling observers come to an end 
within a finite proper time at such singularities. As a con- 
sequence, it is not possible anymore to perform measure- 
ments within GR on such singularities, and an extending 
theory is needed. It is generally believed that quantum 
physics will cure these problems. One usually admits that 
such fundamental theories would possess a classical effec- 
tive theory comprising some of its desired features. In the 
past, candidate theories that soften singularity problems 
were studied [2], but all of these had to modify the vac- 
uum held equations. However, recently Banados and Fer- 
reira discovered that it is possible to avoid the big bang 
singularity by effectively modifying the gravity/matter 
coupling (GraMaC) only [3]. The specific theory pro- 
posed by Banados and Ferreira [3] was coined Eddington 
inspired Born-Infeld (EiBI) gravity. Interestingly, Born- 
Infcld electrodynamics [4] arises as a low-energy effec- 
tive theory of certain string theories [5] , which is one of 
the original motivations to consider Born-Infeld inspired 
structures. 

Our motivations originate somewhere else. Currently, 
the most stringent experimental tests of gravity probe 
spacetime in vacuum but not within dense matter. As a 
consequence, the vacuum field equations are highly con- 
strained, whereas the GraMaC is not. Therefore, there is 
a strong need to better understand implications of such 
modified coupling, both theoretically and experimentally, 
even if the latter is highly challenging. EiBI theory pro- 
vides an optimal playground to understand implications 
of such a modified coupling. The reason is that EiBI cou- 
pled to a perfect fluid exactly corresponds to GR with a 



modified fluid equation of state (EOS). This allows an 
intuitive and explicit understanding of coupling modifi- 
cations, avoiding technical difficulties. This leads to in- 
tuitive new interpretations of the recent results proposed 
in [3, 6]. Implications of modified coupling on energy 
conditions arc illustrated within our framework. More 
important, we reveal the singularity-avoiding mechanism 
discovered in [3] and conjecture a similar effect in more 
general contexts. We further stress that a bimetric-like 
action for EiBI theory provides a clean basis to indepen- 
dently modify the GraMaC and the vacuum field equa- 
tions. 

EiBI as GR with apparent stress. - In this section 
we review the basics of EiBI theory and explicitly show 
the coupling modification in the case of a perfect fluid 
source. EiBI theory introduced in [3] is based on a Pala- 
tini formulation of the following action 
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where g ao is the metric, g = det(g a b), R ao is the Ricci 
tensor built from an independent connection T c ab , and 
Sm is the matter action; in [3], a = 1. Note that most 
previous work on Born-Infeld-like gravity uses a metric 
formulation, see, e.g., [2, 7, 8]. Variation of (1) leads 
to [3, 9] 
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where q ab is an auxiliary metric, q = det(q a b), and 7 = 
a/A. We have absorbed a factor of A in q ao compared 
to [3], both metrics are then identical if T ab = 0. 
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Equations (2) may be written as 

q ac g cb = XS\ - kR\, q ac g cb = r(S a b - 7 kT° 6 ), (3) 

where q ab , g ab are the matrix inverse of q ab , g a b, respec- 
tively and R a b = q ac R cbl T a b = T ac g cb . Combining both 
equations leads to 

R\ = 1tT\ + ^—^S a b , R = 1 rT + 4^-^, (4) 

where R = R a a and T = T a a . The Einstein tensor for 
q ab then follows immediately 

G\[q cd ] = R\ - ^5\R = 7 T a 6 - A<5 a 6 , (5) 

where we defined the apparent stress tensor T a b (the ad- 
dition "apparent" will become clear in the next section) , 

T\ = tT\ + VS\, 1 kV = t-1- ^t 7K T, (6) 

and the cosmological constant A = (A — 1)/k. V plays 
the role of an isotropic pressure addition. Note that r 
can be obtained from T a b by 



T=[dct{8 a b -K 1 T a b )\ 



(7) 



The Bianchi identities G 



ab 



imply T' 



ab 



0, where 



denotes the covariant derivative based on The study 
of EiBI theory now reduces to a study of this potentially 
exotic source TV 

These expressions may still implicitly depend on g ab 
through T a b . However, in important special cases it is 
possible to completely decouple g ab from (5). For exam- 
ple, in the quite general case of a perfect fluid, 

T a b = (p+P) u a u b + PS a b , u a : u a u b g ab = -1, (8) 

it is possible to reexpress T a b in the form of a perfect 
fluid, but in terms of q ab , an apparent fluid velocity v a , 
v a v b q ab = — 1, and an apparent pressure P q and energy 
density p q . In the following, we will refer to g ab , p, and 
P as the real metric, pressure, and energy density, while 
q ab will be referred to as the apparent metric. After 
some algebra using the field equations (2), one finds that 
v a v b = u a u b , while r is computed in a frame adapted to 
the fluid velocity, 

t= [(1+ 7 «P)(1-7^) 3 P- ( 9 ) 
It immediately follows that 

P q = rP + V, Pq + P q = r(p + P), (10) 

and further 

p q =Tp- V, jkV = t - 1 - )pnr^P - p). (11) 

In the case of very dense matter, e.g., within neutron 
stars, the EOS can be approximated to be barotropic, 



p = p(P). Then (10) encodes the apparent EOS. For 
example, dust P — leads to 



3 7 K 
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which smoothly interpolates between adiabatic index 2 
and 1 for small and large, respectively, densities. This 
provides intuitive interpretations for the existence of rel- 
ativistic pressureless stars and singularity avoidance for 
the collapse of dust in the Newtonian limit of EiBI the- 
ory [6]. (The vacuum energy EOS, p = —P, is invari- 
ant.) Yet if a more realistic microphysics is desired, one 
has to find apparent counterparts for other thermody- 
namic variables as well. Fortunately, a scalar quantity n 
(e.g., baryon number density) with associated conserved 
current j a = nu a transforms in an elegant way. Indeed, 
the conservation equation d a \f—gj a = can be rewrit- 



ten as d a \/ '—qn q v a = 0, where n q 



-nrjiv u b ). The 



term —v b u b = yr(l + "fKp) is the time shift between g ab 
and q ab , so the total factor —T/(v b u b ) is the ratio of the 
comoving- frame 3- volume measured with g ab or q ab . 

A minimally coupled scalar field with Lagrangian den- 
sity C m = [(d a 4>)(d a (j))/2 — V] provides a specific micro- 
physics for a perfect fluid. The resulting stress tensor [10] 
admits P = p = + V, where (f> can be thought 
as a parameterizing function. One can then immediately 
compute the equations of motion in terms of P q and p q . 
It should be stressed that the scalar field pressure and 
density implicitly depend on g ab , but it is possible to 
completely eliminate it although leading to very cumber- 
some expressions. 

Observational consequences. Provided one relies 
on gravitational/geometrical imprints measured outside 
matter, one needs to make assumptions on the GraMaC 
to draw conclusions on the source. Remember that EiBI 
theory agrees with GR in vacuum and only modifies the 
GraMaC. If the GR coupling is falsely assumed, then the 
apparent quantities, e.g., the matter EOS, are interpreted 
as real. This is the reason why we qualify these quantities 
as "apparent." Thus by observations of mass-radius rela- 
tions or quasinormal mode frequencies of neutron stars, 
gravitational waves in general, and so on, one may only 
gather information about the apparent EOS. For simi- 
lar reasons, we identify 7 = 8irG here, with the Newton 
constant G. 

On the other hand, the real EOS of very dense matter 
is not fully under control yet [11]. Additionally, probing 
spacctime within such highly interacting matter seems 
to be extremely difficult. However, an important exper- 
imental signature of EiBI theory might be given by the 
case when gravitational waves favor a specific EOS, while 
observations of electromagnetic fields or neutrinos favor 
a different EOS (multi- messenger astronomy [12]). Re- 
member that the latter couple to the real metric within 
the action. Finally, it seems inevitable to constrain the 
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real EOS by experiments on earth (based on electromag- 
netic interactions) and/or by reaching consensus in the- 
oretical predictions. 

It is illustrative to constrain EiBI theory using the 
sun [13], though better constraints are provided in [14]. 
Yet EiBI theory possesses excellent experimental via- 
bility. Notice that the dynamics of compact binaries 
is almost identical to the GR apparent rela- 

tive position and masses exactly equate the real ones. 
Merely subleading EOS-dependent effects arise, which 
appear in post-Newtonian frameworks as, e.g., tidal in- 
teractions [15]. 

Energy conditions. The most important implica- 

tions of matter properties in GR are encoded in various 
energy conditions. These have immediate consequences 
on, e.g., singularity theorems [1], entropy bounds [16], 
and instabilities or other pathologies [17]. The energy 
conditions for perfect fluids are expressed simply by the 
following requirements on pressure and density: 

Null Energy Condition: p + P > 0, 

Weak Energy Condition: p + P > 0, p > 0, 

Strong Energy Condition: p + P > 0, p + 3P > 0, 

Dominant Energy Condition: p > \P\, 

Causal Energy Condition: \p\ > \P\. 

Since T a b is again a perfect fluid, we can check if 
the apparent EOS still fulfills a specific energy condition 
when the real EOS does. In the following, we refer to 
the energy condition expressed in terms of the apparent 
quantities as "apparent energy conditions" and, corre- 
spondingly, those expressed with the real quantities as 
"real energy conditions." 

The apparent Null Energy Condition (NEC) under the 
assumption of real NEC for perfect fluids follows from 
(10) and t > 0. Apart from NEC, energy conditions do 
not generically hold. Indeed, assuming real Strong En- 
ergy Condition (SEC) does not lead to apparent SEC. 
We illustrate this for n > and n < in Fig. 1. As 
can be seen, for k > there is indeed a dark gray region 
where assuming real SEC, we can find apparent SEC vi- 
olation. The SEC ensures by virtue of the Raychaudhuri 
equation that gravity is attractive (for A = 0) and is 
therefore important for certain singularity theorems [1]. 
Yet, singularity avoidance occurs for both signs of k, so 
it is not due to apparent SEC violation. 

One can easily check that Weak Energy Condition 
(WEC) is transported to the apparent sector for posi- 
tive values of k but can be violated for k < 0. Dominant 
and Causal Energy Conditions (DEC/CEC) can be vi- 
olated for both signs of k. This shows the possibility 
of apparent causality violations via superluminal energy 
flows inside matter [1, 16]. (But superluminal speed of 
sound is not excluded even by DEC or CEC.) 

Cosmological case revised. An instructive situa- 

tion considered in [3] is the cosmological case. Depend- 



Apparent SEC violation for real SEC, jk > Apparent SEC violation for real SEC, y*<0 




FIG. 1. Left: Region excluded by real SEC or imaginary 
r (light gray) and where real SEC leads to apparent SEC 
violation (dark gray) for k > in the p — P plane. Right: 
Same plot for k < 0. 
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FIG. 2. Various quantities as a function of time t for k > 
(left) and k < (right) for a dust universe, P — 0. 



ing on the sign of k, the early time, radiation dominated 
universe either freezes or experiences a bounce, while 
the (real) density and pressure remain finite. We have 
redcrived these results. Interestingly, although the real 
quantities are regular at all time, we find that the appar- 
ent sector does have singularities, as predicted by singu- 
larity theorems. 

However, the case of a dust universe P = reveals that 
singularity avoidance is not a generic feature of the the- 
ory. As illustrated in Fig. 2, for k > 0, the Kretschmann 
scalar K signals a singularity in both real and apparent 
sectors, while for n < a bounce with regular real sector 
occurs (analogous to the radiation universe [3]). This is 
opposite of what one might expect from the Newtonian 
collapse [6]. An intriguing explanation is given in the 
following. 

Singularity avoidance. Interestingly, r may be- 

come infinite for finite values of p and P. In the present 
section we argue that this triggers singularity avoidance. 

The first case where r — > oo is given by ^kP — > 1. 
Then both P q /r and p q /r approach a finite constant, 
so P q and p q become infinite. This means that infinite 
P q and pq correspond to finite p and P. More precisely, 
P = I/7K is the maximal real pressure (corresponding 
to a singularity of p q ). Further, the apparent EOS ap- 
proaches P q p q , independent of the real EOS, k, or 
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other thermodynamic variables. A neutron star EOS is 
therefore considerably softened, which qualitatively ex- 
plains the modified mass-radius relation in [6] for k > 0. 
Note that the time shift between apparent and real cos- 
mological times —v b Ub approaches oo. Thus, the appar- 
ent singularity does not form within finite real time. 

Similarly, the second case where r — > oo is given by 
7K/J — > — 1. Again P q /r and p q jr approach a constant, 
but now the apparent EOS approaches P q w — p 9 /3 and 
the maximal (ft < 0) real energy density is p = — I/7K. 
Note that now the time shift goes to zero. 

Whenever one of these two cases applies, gravity is not 
able to produce singular p and P. But it is easy to give 
examples where both cases are avoided, e.g., the dust 
universe with k > explored in the last section. How- 
ever, most realistic EOS will imply that either ^nP — > 1 
(for ft > 0) or 'fK.p — > — 1 (for n < 0) in the high den- 
sity regime. Thus, even when the apparent EOS fulfills 
the SEC and singularity theorems predict a singularity 
in the apparent theory, p and P remain finite in realistic 
situations. This suggests that a regularization analogous 
to the cosmological case might occur for black hole for- 
mation, at least within matter. But further analysis is 
needed to understand under which circumstances world- 
lines of free falling observer will enter the matter and not 
end within finite real proper time, thus avoiding singu- 
larities. Notice that the diverging time shift for k > 
indicates that not even the singularity in the apparent 
sector is encountered by an observer within finite real 
proper time. 

Ghosts, instabilities, and truncations. Recently, 
unitarity analysis was performed for curvature-truncated 
Born-Infeld gravity actions [7]. The result is a generic 
appearance of ghosts in the theory around a constant 
curvature background. These results were found within 
a metric approach, whereas here we follow a Palatini ap- 
proach. Roughly speaking, ghosts are due to a "wrong 
sign" in the kinetic part of the action, which in partic- 
ular is associated with instabilities of the theory. The 
gravitational part of (1) can be equivalently written as a 
linear bimetric-like action, 



So - i 

7 



d A x, 



-(I 



R[q}-2- + -(q ab g ab -2T) 

K ft 



(13) 



which is basically the action in [18] without dynamics for 
g a b. This action is manifestly free of ghosts. 

NEC violation may lead to other types of instabilities 
or pathologies [17]. Assuming real NEC in the full non- 
truncated theory preserves apparent NEC, avoiding this 
type of instabilities. Instead, at the truncated level, it is 
very easy to violate NEC even at second order in trun- 
cation, 

7K 



Pq + P q « (p + P) 



1 



(3P-p) +0[p,P}\ (14) 



p> (2+37K£ 2 )/(47k). One can equivalently truncate the 
curvature sector in the original action (1). This makes 
apparent that truncating a Born-Infeld gravity theory 
can lead to problems, whereas the nontruncatcd version 
is well behaved. 

Truncations also ruin singularity avoidance, e.g., as for 
a post-Newtonian approximation of the apparent source, 



P g =P+ 1 r p 2 + ..., Pq ^p-^-p 2 



(15) 



A straightforward example of NEC violation within the 
truncated theory is P = —p + e 2 , e e K, ft > 0, for any 



Still this recovers the repulsive effect of EiBI theory for 
ft > found in [3, 6]. Indeed, for ft > the apparent the- 
ory possesses both an increased pressure and a decreased 
density compared to the real theory. [But note that the 
K-depcndcnt term in (7) of [3] is coming from the relation 
between q ab and g ab ] 

Conclusions and outlook. - EiBI theory provides a 
prototype of GraMaC modification at high densities in a 
way that avoids the big bang singularity [3] . EiBI theory 
reproduces usual GR in vacuum and to lowest order in n. 
Let us stress further remarkable features of EiBI theory 
discovered in the present Letter and how they translate 
to more general modifications of GraMaC. 

We found that in the case of a perfect fluid, EiBI theory 
exactly reduces to GR, but in terms of an apparent metric 
coupled to a fluid with an apparently modified EOS. The 
notion of apparent quantities generalizes to more generic 
GraMaC, although real fluids will not always map to ap- 
parent fluids. On the one hand, this presents a challenge 
for experimental constraints on the theory (and concur- 
rently on the real nuclear EOS from neutron star ob- 
servations), since it is very difficult to probe spacetime 
within matter. On the other hand, this demonstrates 
the excellent current and future experimental viability of 
modified GraMaC theories. Notice that EiBI theory is 
viable also from a theoretical point of view, as it appears 
to be free of instabilities, ghosts, or other inconsistencies. 
Such theories can also make regular matter appear exotic 
as illustrated by the energy conditions. 

The similarity to GR allows for a transfer of the vast 
amount of knowledge on GR to EiBI theory in the future, 
e.g., usage of existing codes for full numerical simulations, 
exploration of the Oppenheimer-Snyder solution [19] for 
the collapse of dust (where now P q = 0), or study of 
naked apparent singularities as hypothetical astrophys- 
ical objects [20]. This should also be the case in more 
general situations since essentially the description of the 
source needs to be rcimplcmcntcd. 

We revealed that singularity avoidance in EiBI the- 
ory occurs when a pole in r, Eq. (7), is approached and 
illustrated this mechanism by singular and nonsingular 
examples in the cosmological case. Though in the nonsin- 
gular case the apparent curvature becomes singular, this 
is completely harmless for measurable (real) quantities. 
It is most interesting to study explicitly the geodesies 
or the validity of entropy bounds during a collapse to a 
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black hole in the future. It can be expected in some par- 
ticular cases EiBI that allows physical predictions where 
GR fails. 

An intriguing new interpretation of the gravitational 
interaction derives from (13), in particular when viewed 
as a rather unusual effective action. This action consists 
of the "bare" Einstcin-Hilbcrt part for the dynamical 
metric q a b, which induces a dynamics to the measurable 
(matter-coupled) metric g a b through a constraint-type re- 
lation. The "cutoff" (7ft) -1 enters (13) as a coupling con- 
stant. Instead, in bimetric theories, usually both metrics 
are truly dynamic. This makes the action (13) interest- 
ing as a basis for independent modification of GraMaC 
and vacuum field equations, opening up a new class of 
modified theories of gravity. 
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